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Abstract. The mass- loss properties of carbon AGB stars are not very well 
constrained at present. A variety of empirical or theoretical formulae with differ- 
ent parameterisations are available in the literature and the agreement between 
them is anything but good. These simple mass-loss prescriptions are nonethe- 
less used in many models of stellar evolution without much consideration of their 
applicability in various cases. We present here an on-going project aiming at a 
better description of the mass loss, that could be used to improve stellar evolu- 
tion models - especially the evolution during the TP- AGB phase. As a first step, 
we have considered the mass-loss evolution during a He-shell flash. Using stellar 
parameters sampled from a stellar evolutionary track, we have computed the 
time evolution of the atmospheric layers and wind acceleration region during a 
flash event with detailed frequency-dependent radiation-hydrodynamical models 
including dust formation. We find that existing simple mass-loss prescriptions 
imply mass-loss evolutions different than our model. Based on these results, we 
have also simulated the subsequent long-term dynamical evolution of the circum- 
stellar envelope (CSE), including the formation of a detached shell. The second 
step of the project deals with the dependence of mass loss on the basic stellar 
parameters. At the moment we are computing a large grid of wind models for 
C-rich AGB stars. Preliminary results show that simple parameterisations are 
difficult to obtain in certain regions of the parameter space considered, due to 
strong non-linearities in the wind mechanism. 



1. Introduction 

Models of mass loss on the AGB are essential for understanding the late stages 
of stellar evolution. In practice, the modelling of stellar evolution requires a 
prescription of mass-loss rates as a function of basic stellar parameters. Empir- 
ical studies, such as the classical one by Vassiliadis & Wood (1993), can only 
provide a limited amount of information, since one cannot easily study the effect 
of individual parameters. The reason for this is of course that the evolution of 
the stellar parameters is interconnected and we can only observe stars at a given 
evolutionary stage. Observations can tell us about the mass loss for a given 
set of stellar parameters, but the exact dependences on these parameters are 
unfortunately degenerated. 

Numerical modelling has one major advantage here: it can be used as a tool 
for constraining the actual physics involved in wind formation. From the early 
efforts, e.g. by Bowen (1988) and Fleisch er et al. (1992), to our latest models 
( Hofner et al. 120031 ; iMattsson et al.l [20071 ) . a lot of the physics behind dust- 
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Figure 1. Formation and dynamical evolution of a detached shell: snap- 
shots of the velocity and density structure according to the CSE model. Upper 
panels: the velocity field; lower panels: the logarithmic density. Left: piston 
amplitude Av p — 4.0 km s _1 ; right: piston amplitude Av p — 6.0 km s _1 . 



driven winds has been brought to light. Much work still remains to be done, but 
we have now reached a reasonable level of realism in the modelling, as shown 
by detailed comparison with observations. A picture has emerged with dust- 
driven winds as a complex phenomenon, where the onset of mass loss requires 
particular circumstances. In fact, C stars with strong dust-driven winds may be 
a rare species even if most C stars apparenly lose mass. With this in mind, one 
may ask: What are the consequences for stellar evolution, nucleosynthesis, and 
consequently the chemical evolution of galaxies? 



2. Mass-Loss Modulations during the TP-AGB Phase 

The common hypothesis about how detached shells around TP-AGB stars are 
formed is that these structures are a consequence of the mass-loss modu lations 
during the phase of He-shell flashes/thermal pulses ( Qlofsson et al.|[l990l ). Pre- 



vious modelling has shown that interactions between different wind phases that 
arise from variatio ns in the wind velocity and m ass-loss rate are capable of pro- 
ducing such shells ( Steffen Sz Schonbernerl [20001 ) . 



Using the RHD code for dynamic atmospheres by Hofner et al. (2003) in 
combination with a stellar evolu tionary track of a M = 2 M Q , Z = 0.01 star 
as input (see iHerwid |2000| . 120041 . for details about the stellar e volution model) 



we h ave calculated the wind evolution during a He-shell flash (jMattsson et al 
2007). With this evolution computed, we constructed a time-dependent inner 
boundary condition for a hydrodynamic model of the circumstellar envelope 
(CSE) evolution. Th e observed wind and mass-loss properties associated with 



(^SHjJ evolution, ine observed wind ana mass-loss prop erties associated witn 
detached shells (e.g. lOlofsson et aD l2000l : ISchoier et ail 120051 ) are rather well 



reproduced using stellar parameters sampled from the evolutionary track. We 
tried different amplitudes for the piston boundary condition which is used to 
simulate the stellar pulsation and found that too high amplitudes lead to very 
little wind interaction and thus no shell structure is formed (see Figure [1]) . Too 
low amplitudes, on the other hand, tend to prohibit any wind formation at all. 
It seems that for the stellar parameters investigated here, the optimal piston 
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amplitude is ~ 4 km s _1 , and this value may perhaps be used to constrain the 
pulsation properties. 

It is also interesting to note that existing mass-loss prescriptions for AGB 
stars imply mass-loss evolutions that are not only different in compar ison with 



our d etailed model, but are rarely consistent with each other, either (see lMattsson et al 



20071 . for further discussion). Furthermore, constra i nts on the mass-los s evolu- 



tion derived from observations ( Olofsson et alJl200Cl : ISchoier et al.ll2005l ) appear 
to agree better with our results than with these mass-loss formulae. 



3. Mass Loss as a Function of Stellar Parameters 

We are currently in the process of computing a large grid of dynamical C star 
atmosphere models. The overall purpose of the grid is to explore how the mass- 
loss rate depends on the stellar parameters. In a future step of this project we 
will consider how pulsations and winds affect the spectral properties of C stars 
as well, but here we will discuss only the wind formation and the mass loss itself. 

Before summarising the numerical results, it is worthwhile to consider some 
fundamental aspects of the onset of dust-driven winds from a basic physics point 
of view. The condensation of dust grains is very sensitive to temperature and 
is prohibited by high gas temperatures and low densities, which means that 
every C star has a condensation radius (R c ) - a critical radius inside which dust 
grains cannot condense - as well as an outer radius where the condensation 
stops. Furthermore, in order to levitate the material out to R c , the acceleration 
of the gas due to pulsations and radiation must be large enough. There cannot 
be any circumstellar dust formation and hence no dust-driven wind if insufficient 
material is transported out to R c . The location of R c will depend on the effective 
temperature of the star, as well as on its photospheric radius. Then, quite 
obviously, there must be a temperture-dependent threshold for the formation 
dust-driven winds. The sudden drop in the mass-loss rate and wind velocity 
above T e g ~ 3000 K is probably even steeper than Figure [5] seems to indicate, 
due to a grid spacing of 200 K. 

Another very obvious threshold emerges as we consider the wind accel- 
eration process after the onset of dust formation. Beyond the condensation 
radius R c , the criterion for wind formation is that the radiative acceleration 
must be strong enough to overcome the gravitational field of the star, i.e. 
Tdust > r cr jt, where T is the ratio of the radiative to the gravitational accelera- 
tion (see e.g. lDominik et aL 1990). The radiative acceleration is proportional to 



the dust opacity Kdust, which is a function of the abundance of free ca rbon eq 
and the mean degree of dust condensation f c (cf . iHofner &: Dorfil Il997l ) . Since 
^dust oc we have that 

Pgas 

r du st oc — , = ■ ec/c = — 7 ljeo/c. (1) 

Thus, for a given mass and luminosity, there must exist a critical value of C/O 
(or more exactly eq — £o) fo r which T dust = T cr i t . We also note that this critical 
C/O cannot be arbitrarily close to unity for realistic stellar parameters. 

Let us now turn to the numerical modelling results. We present here a sub- 
grid of models of stars with L± = 10 4 L , M = 1 M , Z = Z & and Av p = 4.0 
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km s _1 , where we vary T e g and C/0. To eliminate the pulsati on period (free pa - 
rameter), we employ an empirical period-luminosity relation ([Feast et al.l fl989). 
The predicted threshold appears as expected and we find that for C/O < 1.2 
and/or T e g- > 3200 K no dust-driven wind is formed with the other parameters 
fixed as above (see Fig. [2]). Thus, any mass loss in this part of parameter space 
is likely due to some other mechanism. Since the surface gravity scales with 
T e g, it is quite expected that mass loss scales with T e g, too. But we also see a 
rather strong dependence on C/O for both the wind velocity and the mass- loss 
rate, which is quite interesting in comparison with previous studies of this kind. 
Arndt et al. (1997) as well as Wachter et al. (2002) argue for a weak dependence 
on C/O, which stands in sharp contrast to the results presented here. This is 
probably due to the parameter range for which these formulae where developed, 
as well as differences in the basic model assumptions. 





Figure 2. Mass loss in M Q yr 1 (left) and wind velocity (right) as functions 
of T eff and C/O for a C star with L„ = 10 4 L , M = 1 M Q , Z = Z and 
Aw p = 4.0 km s" 1 , according to our detailed numerical wind model. 



4. Effects on Stellar Evolution and Nucleosynthesis 

Mass loss affects stellar evolution is various ways. For C stars the duration of 
the AGB phase and the number of thermal pulses is almost solely determined by 
the mass-loss rate. The internal structure depends on the mass-loss rate as well 
( Blocker! [l995). which in turn affects the fundamental stellar parameters. Since 
the mass-loss rate depends on these stellar parameters, there is an interesting 
feedback mechanism at work here, which means that the mass-loss prescription 
put into a stellar evolution model is critical. The observed initial-final mass 
relation is our best empirical constraint, but it can just tell us the total mass 
lost during the RGB and AGB phase. Hence, much of the AGB evolution will 
only be guess-work as long as we do not understand mass loss properly. 

Perhaps even more interesting, in the context of this conference, is the effect 
on nucleosynthesis. It is well known that varying the mass- loss rate can have pro- 
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found effects on the chemical yields for AGB stars (see e.g. Ivan den Hoek &: Groenewegen 
19971 ). As we pointed out above, the mass-loss rate determines the number of 



thermal pulses, which represents the maximum possible number of third dredge- 
up (TDU) episod es. Furthermore, it can also determine if and when the TDU 
begins and ends (jKarakas et al. 2002). Needless to say, a correct implementa- 
tion of mass loss is essential to obtain realistic yields for AGB stars, and the 
AGB contribution to the chemical evolution of galaxies therefore still has to be 
regarded as uncertai n. This also become s evident when different sets of AGB 
yields are compared ( Qlofsson et al.|[2007l ). 



5. Concluding Remarks 

Existing mass-loss prescriptions for AGB stars are inconsistent with each other 
and simple mass-loss prescriptions might be hard to obtain. Here we have tried 
to show that it may be dangerous to use parametric mass-loss formualae includ- 
ing too few stellar parameters. Moreover, the first results from our new detailed 
model grid suggests that C stars with strong winds may be quite rare. How 
would this affect stellar evolution? Nucleosynthesis? Chemical evolution? The 
present dynamic atmosphere models compare well with observations and it is 
therefore likely that we are beginning to understand the wind mechanism well 
enough to be able to reduce the current uncertainties related to the mass-loss 
evolution of C stars. 
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Discussion 

Speck: You state that for C/O < 1.2 mass-loss doesn't happen. But the literature 
suggests that C/O < 1.2 for most C stars. And most C stars have dust (SiC 
11 ixva feature). How do you reconcile these? 

Mattsson: Yes, about 75% of all observed C stars have C/O ratios below 1.2, 
but still mass-loss is detected. What our models tell us is that below 1.2 there 
is no strong, dust-driven wind. The observed mass-loss must be due to some 
other mechanism, if the dust is there. 

Gallino: Maybe combining with > 1.5 M Q models is a better choice. For 1 M 
no TDU is typically found. 

Mattsson: Yes, this is just a first result. The grid will be vastly extended to 
cover all relevant parameter ranges. 

Willson: You slid by the piston amplitude issue. Good piston constraints are 
essential for a good theoretical M relation. We found the following very impor- 
tant: Keep track of the work done on it by the piston. Limit this to be < L. 
See Cox & Giuli. This gave d log M/d log L consistent with Vassiliadis & Wood, 
the only reliable empirical determination. 

Mattsson: Very well. If the pulsation were too strong we would see this in the 
luminosity variations, but we don't. 

Aringer: I also see a problem that there is a large number of carbon stars with 
higher temperature and lower C/O losing mass. Maybe dust is not the driver 
there. However, do not trust any published C/O ratios for AGB variables. There 
are large uncertainties. 



